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A  METHOD  TO  DETERMINE  MOISTURE  PROFILES  FROM  TOTAL 
MOISTURE  WEIGHT-GAIN  DATA  IN  POLYMERIC  COMPOSITES 


By 


Y.  Jack  Weiisman 

The  University  of  Tennessee  at  Knoxville 


Abstract 


This  article  presents  a  numerical  approximation  scheme,  based  upon  the 
method  of  moments,  to  evaluate  moisture  distributions  within  a  com¬ 
posite  plate  from  data  on  total  moisture  weight-gain.  Although  it  is  as¬ 
sumed  that  the  diffusion  process  remains  linear,  the  scheme  applies  to 
moisture  uptake  data  which  depart  from  classical  predictions.  Two 
examples  are  provided  to  demonstrate  the  validity  of  the  scheme. 


1.  INTRODUCTION 

The  most  commonly  utilized  technique  to  record  moisture  uptake  in 
polymeric  composites,  and  materials  in  general,  is  to  expose  them  to  an 
ambient  environment  and  record  weight  gain  with  time.  However,  for 
most  engineering  purposes  it  is  important  to  know  the  spatial  dis¬ 
tributions  of  moisture  rather  than  the  integrated  values  which  corre¬ 
spond  to  weight  gains.  Several  experimental  techniques  were  employed 
to  record  moisture  distributions^'^^^^  but  these  are  rather  cumbersome, 
to  say  the  least. 

When  moisture  weight-gain  data  under  exposure  to  constant  ambient 
conditions  follow  the  values  predicted  by  classical  diffusion  theory,  as 
shown  in  Fig.  1,  there  is  good  reason  to  believe  that  the  moisture 
distributions  can  be  "back  tracked"  in  accordance  with  well  known 
classical  results^-^  It  may  be  noted  that  the  computations  of  these  dis¬ 
tributions  are  somewhat  tedious  and  that  the  scheme  proposed  in  the 
present  article  appears  to  be  more  efficient.  Nevertheless,  for  M(t) 
which  agrees  with  Fig.  1,  the  determination  of  spatial  distributions 
follow  well  established  expressions. 

In  many  circumstances  weight  gain  data  deviate  rather  significantly 
from  the  classical  predictions  shown  in  Fig.  1.  These  departures  w-ere 
studied  extensively,  mainly  by  polymer  scientists,  and  attributed  to  a 
variety  of  causes  and  mechanisms.  One  plausible  interpretation  at¬ 
tributed  those  departures  to  the  time  dependent  configurational  rear¬ 
rangements  w'ithin  the  polymer  chains,  akin  to  viscoelastic  mechanical 
response  (see  ref.  (4)  for  a  listing  of  some  of  these  works).  According  to 
this  interpretation,  the  chemical  potentials  of  the  polymer  and  the 
surroundings  do  not  equilibrate  instantaneously  even  under  constant 
ambient  conditions.  The  gradual  approach  to  equilibrium  can  be  ex¬ 
pressed  by  a  time-dependent  boundary  condition  within  the  context  of 
an  otherwise  linear  ("classical")  diffusion  boundary  value  problem. 
Since  the  boundary  condition  is  generally  unknown,  it  is  impossible  to 
"back  track"  moisture  distributions  from  weight-gain  data  using  text¬ 
book  solutions.  The  scheme  developed  in  this  article  accomplishes  this 
task.  However,  it  should  be  noted  that  the  precision  of  the  predicted 
moisture  distributions  cannot  exceed  accuracy  of  the  weight-gain  data. 

2.  ANALYSIS 

Consider  a  "classical,"  one  dimensional  diffusion  process  which  is  gov¬ 
erned  by  the  field  equation 
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t  >  o. 


-L  <  X  <  L 


(1) 


at 


=  D 


a  c 


dx‘ 


where  c  is  the  concentration  of  the  solvent  (water)  in,  say,  percent 
weight-gain,  D  the  diffusion  coefficient  in  m^/sec.,  t  denotes  time  and  x 
is  distance  across  a  plate  of  thickness  2L.  Assuming  an  initially  dry 
plate,  we  have 


c(x,o)  =  0  -L  <  X  <  L. 
The  boundary  conditions  are  taken  to  be 

c(±L,t)  =  f(t) 


(2) 


(3) 


where  f(t)  reflects  the  time-dependent  response  of  the  polymeric 
component  within  the  plate.  This  time-dependence  is  assumed  to  occur 
even  under  constant  ambient  conditions,  accounting  for  total  weight 
gains  that  depart  from  predictions  based  upon  f(t)  =  constant. 


While  most  experimental  data  give  total  weight  gain,  namely 
M(t)  =  I  c(x.t)  dx  the  detailed  distribution.  c(x,t).  is  of  higher  rele¬ 
vance  for  computations  of  residual  hygro-thermal  stresses  in  composite 
laminates.  The  following  approximate  scheme  enables  us  to  unravel  the 
information  contained  in  M(t)  to  obtain  c(x,t).  The  approach  extends  the 
"method  of  moments"  developed  for  unbounded  regions^-*'^ 

The  formulation  presented  in  equations  (l)-(3)  gives  c(x,t)  =  c(-x,t).  | 
This  evenness  in  x  suggests  the  construction  of  the  even  mo  nents  de¬ 
fined  by 


MCpXt) 


c(x,t)  dx 


For  future  use  define  the  k-fold  time  integral  of  M(t) 


(4) 
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JpTk=l  rtk-l  /T2  f1\ 

I  . I  I  M(i)  dx  dxi  . dxk-i  (5) 

0  Jo  Jo  Jo 

and  the  k-fold  time  integral  of  boundary  condition  multiplied  by  2L 

Jptk=l  /tk-1  ^t2  pi 

I  .  I  I  c(L,t)  dx  dX]  . dxk-]  (6) 

0  Jo  Jo  Jo 

Note  that  M(0)(t)  =  l(0)(t)  =  M(t)  and  J(0)(t)  =  2LC(L,t).  Furthermore,  for 

c(±L,t)  =  Cl  =  constant  we  have  2Lcl  =  M„  whereby,  for  time-indepen¬ 
dent  boundary  condition,  we  have 

J(k)(t)  =  M,.  tVk! 

The  method  of  moments  derives  from  the  fact  that  integration  of 
equation  (1)  yields 
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=  2  DL^P  -  2p  2  DL^P-’  c(L,t) 

dx  x=L 

+  2p  (2p-l)  D  j^j^x2<P->)  c(x,t)  dx 

=  L2Pl(0)(t)  -  2pDL“<P-*^  J(0)(t)  +  2p(2p-l)DM2(p-i)(t)  (8) 

Integration  with  respect  to  t  and  employment  of  mathematical  induc¬ 
tion  yield 


M(2p)  (t)  =  l'p  1(0)  (t)  +  I 

^  ^  k=i  (2p-2k)! 


(2p)!_  L2(p-k)  j 


J(k)(t)  I 
2(p-k)+lJ 


Note  that  equation  (9)  involves  the  weight-gain  data  M(t)  and  their  time 
integrals,  a  procedure  which  is  much  less  susceptible  to  errors  than  data 
differentiation. 

.At  this  stage  it  is  expedient  to  introduce  the  non-dimensional  coordinate 
c  =  x/L,  with  a  corresponding  distribution  c  (q,t) .  Consequently,  we  have 


M(:p)(t)  = 


=  1  ^'Pc(;,t)d^  =  M(2p)(t)  /  L“P^’ 


(10) 


Expand  c(c,t}in  a  series  of  even-order  Legendre  polynomials  P2j  (4) 


c  (tt)  =  X  a2j(t)H(t)P2j  (^) 
.1=0 


-!<£<! 


(11) 


where  a^jd)  are  the  yet  unknown  coefficients  and  H(t)  is  the  Heaviside 
unit  step  function. 

In  view  of  the  orthogonality  of  the  Legendre  polynomials  we  have^^^ 


c  (^.t)P2k(^)d^  =  a2k(t)H(t) 


4k  +  l 


(12) 
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Substitution  of  equation  (11)  into  equation  (12)  and  employment 

of  equation  (10),  determine  a2k(t)  in  terms  of  the  moments  M(2j)(t) 
j  =  0,1,  ...  k  as  follows 


^  a,,(0  = 

4k+l  (2k-l)!!(2k):! 


M(2kK0 


(4k-3)!! 

(2k-3)!!(2k-2)!!2!! 


M(2k-2KO 


(4k-5)!! 

(2k-5)!!(2k-4)!!4!! 


M(2k-4Kt)  -  ....  +  (-1) 


k  (2k-l)!! 
(2k):! 


(13) 


In  equation  (13)  we  used  the  notation  (2q-l)!!  =  1-3-5'  ....  •(2q-l)  and 
(2q)!!  =  2-4-6  ....  •(2q). 

For  instance 

2a()(t)  =  M((i){t).  -  aMi)  =  ^  ~  M(o)<t). 

“  2  2 

2-  a4(t)  =  .^^3  5M(4){i)  -  30M(2)(t)  +  3M(())(t)], 

9  8 

-2-  af,{i)  =  -2-  23  1M(6)(t)  -  315M(4)(t)  +  105.M(2Ki)  -  5M(oKi)  >  etc. 

13  1  6  ‘ 

In  view  of  equation  (13)  it  is  obvious  that  the  representation  (11)  suf¬ 
fers  from  accumulations  of  truncation  errors.  These  errors  are  due  to 
the  inherent  scatter  in  the  data  for  M(0)(t)  as  well  as  computational  im- 
precisions  in  the  numerical  integrations  included  in  the  higher  order 
moments  M(:k)(t),  k=l,  ...  p.  The  examples  presented  in  the  foregoing 
section  indicate  that  very  good  results  for  c(x,t)  can  be  obtained  by 
taking  k  <  4,  in  which  case  these  truncation  errors  may  not  exceed  a  few' 
percentage  points. 

3.  COMPUTATIONAL  RESULTS 

Computations  were  performed  for  two  test  cases  where,  for  verification 
purposes,  both  c(x,t)  and  M(t)  are  known  analytically.  In  both  cases  we 

took  L  =  1.0  mm  and  D  =  5.8  x  10'“^  mm^/sec.,  while 
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case  (a)  :  c(±L,t)  =  1.0%  (constant) 

case  (b)  :  c(±L,t)  =  0.8  +  0.2  (3  =  l/200. 

Expressions  for  M(t)  and  c(x,t)  are  given  in  Ref.  (7).  The  values  of  M(t) 
were  employed  as  simulated  data  to  generate  M(2p)(t)  according  to 
equation  (9)  and  recreate  c(x,t)  according  to  equation  (11)  by  solving  for 
a2j(t)  in  accordance  with  equation  (13). 

Results  are  shown  in  Figures  2  and  3.  These  figures  show  distributions 
of  c(x,t),  employing  up  to  six  terms  in  the  series  (11),  for  various  times. 
The  numerals  accompanying  the  various  curves  indicate  the  number  of 
terms  in  the  series,  equation  (11).  employed  in  their  evaluation.  For 
comparison  purposes,  the  exact  distributions  c(x,t),  which  are  known  for 
cases  (a)  and  (b),  are  shown  in  solid  lines.  Note  that,  for  the  cases  under 
consideration,  very  good  results  were  obtained  by  using  no  more  than 
three  or  four  terms  in  the  series  (11).  The  divergence  that  is  caused  by 
taking  additional  terms  is  due  to  trucation  errors. 

4.  CONCLUDING  REMARKS 

The  computational  scheme  developed  in  this  work  employs  total 
weight-gain  data  to  determine  moisture  distributions  across  the  thick¬ 
ness  of  polymeric  composite  laminates.  The  scheme  is  most  useful  for 
cases  where  weight-gain  data  depart  from  predictions  of  classical,  linear 
diffusion  theory,  when  no  "textbook"  solutions  are  available  to  relate 
total  uptake  to  spatial  distributions.  It  should  be  noted  that  the  scheme 
is  based  upon  the  premise  that  linear  theory,  as  expressed  in  equations 
(1)  and  (2),  still  applies,  and  that  the  sole  cause  for  deviation  from 
classical  weight-gain  is  the  time-dependent  material  response  which  is 
reflected  in  equation  (3).  In  a  forthcoming  work  we  intend  to  propose 
some  criteria  to  examine  if  moisture  uptake  data  concur  with  the  above 
premise  and  devise  a  method  to  determine  f(t).  In  addition,  it  should 
be  noted  that  the  diffusion  coefficient  D  is  also  generally  unknown  and 
must  be  determined  experimentally.  It  is  our  intention  to  address  this 
aspect  as  well. 
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The  total  weight  gain  ratio  M(t)/M(oo)  vs.ft*  according  to  Pick's  Law.  with 
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locations  where  departures  from  straight-lines  exceed  1%.  t*  =  Dt/L  . 
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FIGURE  2.  Moisture  distribution  profiles  c(x,t)  vs.  distance  x  computed  according  to  equation  (11), 
with  k  indicated  on  each  dashed  curve.  The  exact  solution  is  shown  by  solid  line. 

Cl=1,  L=1,  D  =  5.8  X  10  ^  (a)t  =  50,  (b)t  =  100,  (c)t  =  200,  (d)t  =  1000. 
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FIGURE  2.  Moisture  distribution  profiles  c(x,t)  vs.  distance  x  computed  according  to  equation  (11), 
with  k  indicated  on  each  dashed  curve.  The  exact  solution  is  shown  by  solid  line. 

Cl=1.  L=1,  D  =  5.8  X  10-^  (a)t  =  50,  (b)t  =  100,  (c)t  =  200,  (d)t  =  1000. 


Moisture  Distribution  Profile 


nGURE3:  Same  as  Figure  2,  except  that  Cl=  0.8  +  0.2(1  -  exp(t/200)l 
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FIGURE  3:  Same  as  Figure  2,  except  that  0^=  0.8  +  0.2|  I  -  exp(t/200)I 
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FIGURES:  Same  as  Figure  2.  except  that  Cj  =  0.8  +  0.2[1  -  exp(t/200)) 


